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Abstract 

We study the Aharonov-Bohm effect for an exciton on a nano-ring using a 2D attractive fermionic 
Hubbard model. We extend previous results obtained for a ID ring in which only azimuthal 
motion is considered, to a more general case of 2D annular lattices. In general, we show that the 
existence of the localization effect, increased by the nonlinearity, makes the phenomenon in the 
2D system similar to the ID case. However, the introduction of radial motion introduces extra 
frequencies, different from the original isolated frequency corresponding to the excitonic Aharonov- 
Bohm oscillations. If the circumference of the system becomes large enough, the Aharonov-Bohm 
effect is suppressed. 
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I. INTRODUCTION 



Progress in micro-structuring technology allows the construction of nano-sized InGaAs 
rings by self-assembly l| ^[J or using lithographic techniques j^,!^. Similarly PbSe-based 
nano-rectangles have been synthesized in solution through oriented attachment of nanocrys- 
tal building blocks p. In these systems, electrons and holes can propagate coherently and 
the existence of bound states — so-called excitons — offers an opportunity to explore the 
Aharonov-Bohm (AB) effect Q,!^. 

netic flux due to its finite sized Bohr radius This effect, in the framework of excitons 
in nano-rings, has been studied previously by means of simple models, such as a continuous 



The exciton, although neutral, can be sensitive to a mag- 
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BH HQ, Q, 



or in continuous two-dimensional 



one-dimensional (ID) quantum ring 
(2D) models for which some numerical solutions can be found 114. The effects of an addi- 

n 

tional external electric field have been also been studied [15|. In all cases, it has been shown 
that AB oscillations can exist in a finite-width nano-ring, but when the ring becomes large 
enough, the effect is suppressed. A similar effect has also been predicted for bi-excitons . 

In this paper we consider lattice models of such nano-rings with an electron and hole 
subject to a perpendicular and uniform magnetic flux. Our models allow for a finite number 
of transport channels, either arranged annularly, i.e., coupled rings laid out in a 2D plane 
with increasing radius, or stacked vertically. Therefore, we can study in a controlled way 
how the excitonic AB oscillations evolve when additional rings are added, either in-plane 
or vertically. We describe the electron-hole system on the rings by means of an attractive 
Hubbard model where electron and hole are modelled as spin-up and spin-down particles 
with different effective masses. We focus our interest on the study of the AB oscillations, 
their dependence on the size of the rings, and on different combinations of lattice parameters. 

For attractive Hubbard models, the existence of bound states of electrons and holes - 
identical to spin-singlets when using the standard electronic interpretation — is well known 
16], and has been previously also been studied in the framework of quantum breathers 
17| . In these cases, if the anharmonicity (interaction strength) is strong enough, there is an 
extended Bloch state with two or more particles in a strongly correlated state. There exist 



various theoretical 
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20H and numerical 121 



of these states in different quantum systems (22 , 



results as well as experimental observations 



Our work is organized as follows: in the next section we present the model. In Section 
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IIHI we study the ID model and analytical expressions are obtained and compared to the 
ID continuum approach. In Section HVl we consider a simple plate-shaped 2-ring model and 
compare the AB oscillations with those obtained in the ID case. In Section we extend 
our results to a finite number (> 2) of rings. In Section IVH we consider a model in which 
the rings have been stacked vertically, so that the same flux is threading each ring. Finally, 
in Section fVIIl we summarize our findings and present our conclusions. 



II. THE MODEL 



We consider an anharmonic 2D lattice of N rings, each ring with M sites, subject to 
a uniform magnetic field perpendicular to the lattice, as shown in Fig. ^ On this lattice 
we consider an electron and a hole. We assume that the interaction potential is short- 
ranged and the system can be described by an attractive fermionic Hubbard model. Also, 
we will consider only the possibility of hopping between nearest-neighbour sites. Thus, the 
Hamiltonian is given by 

N M 

H = ~~ 7 ^ ] ^ ] a n,m a n,mbn jm ^n,m~ (1) 



n=l m=l 
,V-J M 

, : m u n+l,m ~T~ a n+l,m a n,m J t~ 



£ £ 4 [<: 



71=1 771=1 

N M 



/ j / j L n I e I u n,m u "n,rn+l T H>U n ,m+l u n,m J T" 

71=1 771=1 

-2ni<p n /M f t „ _i_ w/it h 

c I u, n,m+l u "ri,m < h lu n,m u n,m+l 



here a^ m (a n)m ) and &] vm (& n ,m) are raising (lowering) operators for electrons and holes, re- 
spectively, satisfying the standard fermionic anti-commutation relations. The parameters tn 
represent the hopping coefficient between neighbouring sites along the n th ring, and the 
hopping coefficients between neighbouring sites, one corresponding to the nth ring and the 
other one to the (n + l)th ring (in absence of magnetic field). In general, these coefficients 
will depend on the distance between sites. For an annular ring structure as in Fig. the 
site-to-site distances increase as one moves further away from the centre. In an single ring 
with no interaction, (7 = 0), the effective mass of the electron and hole at the extrema of the 
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FIG. 1: Sketch of the model. The lattice is a 2D system of N rings and M sites in each ring. The 
system is subject to an uniform magnetic field perpendicular to the lattice, with the magnetic 
flux. 

conduction and valence band are proportional to the inverse of the square of the distance 
between sites [12]. In our model, we have chosen, for a fixed value value of the number of 
sites, a fixed effective mass for all the rings with different radius and, therefore, we model 
the intra-ring nearest-neighbour hopping coefficient by assuming a dependence according 
to an inverse-square law. Also, we consider a similar dependence for the inter-ring hopping 
coefficient. Thus, 

n (r„ +1 -r n ) 2 ' » 4r2sin 2 (7r/M)' {) 
where r n is the radius of the nth ring and d n = 2r n sin(7r/M) is the chord distance between 
two neighbouring lattice points of the nth ring. Also, the formal continuum limit can be 
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obtained by taking this distance going to zero. We will suppose that the distance between 
rings is constant, = t x , and thus independent of the ring. The parameter 7/e represents 
the ratios of anharmonicity (interaction parameter) to nearest-neighbour hopping energy, 
\i the ratio of the effective masses of electrons and holes at the bottom and the top of the 
conduction and valence band respectively, and (p n = <& n /$>o, the magnetic flux through 
;he nth ring and $0 the flux quantum. In general we consider e = 1 and u = 0.2, as in Ref. 



3. 



The ID Hubbard model with nearest-neighbour hopping terms has been solved ex actly 
in the S e mi nal papers by L.eb and WnH and with flux hy Shastry and Sutherland HQ- 
Recent results are reviewed in Ref. [28j . Similarly exact results for extensions of the model 
are rare. An approximation to calculate the spectrum of the Hubbard model has been 
proposed, see, e.g. 25], where the exact spectrum and thermodynamics for a long-range 
hopping Hubbard chain with linear dispersion is calculated. In this paper, we will use the 
number-state method to calculate the eigenvalues and eigenvectors of the Hamiltonian 
operator in the sector with a single electron and a single hole. We use a number-state 
basis 
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where n e (m _ 1)N+n {n\ m _ x)N _ 
ring m. We always have 



n ) represents the numbers of electrons (holes) at site n of the 
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1 and $2 
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1, 



1. A general wave 



function is |^/) = ^2^ n y C{ n }|^{ n }). It is possible to block-diagonalize the Hamiltonian op- 
erator using eigenfunctions of the periodic translation (or rotation) operator T defined as 



n,rn+ 



1 T). In each block, the eigenfunctions have a fixed value of 



T6 t = fe f 

n,m n,m+ 

the total momentum K, with r = exp(zi^) being an eigenvalue of T such that K = 2i7r/M 
and £ integer [l7|. 

In general, if the anharmonicity parameter is large enough, there exists some isolated 
eigenvalue for each k which corresponds to a localized eigenfunction. By "localized", we 
mean there is a high probability for finding the two fermions at the same site. This is 



illustrated by the band structure in Fig. which shows two isolated bands below the two 
"continuum" bands. Due to the rotational invariance of the system, there is the same 
probability of finding two fermions at any site of a ring. In these cases, analytical expressions 
can be obtained in some asymptotic limits 




FIG. 2: Band structure of the eigenvalues E(k). M = 25, N = 2, 7 = 4, ip\ = Lp 2 = 0, d\ = 1, 
e = 1, pb = 0.2. Ar = r 2 — r\ = 0.5. 
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III. A SINGLE RING 



A. The Hamiltonian matrix 

In the simplest case, N = 1, the system reduces to a ID ring of M sites and radius r. 
The Hamiltonian matrix in each momentum subspace is given as 

7 q* . . q 
q q* . 
q q* . . 



(K) 



(4) 



. . . q q* 
q* . . . q 

where q = t^e 2m(p ' M (/i + e _iK ) and we have omitted the index labelling the ring number. 
In the case of no interaction (7 = 0) between the electron and hole, the energy of the system 
is given by 



E. 



k e ,k h 



-2t 11 cos ( k e 



2tt$ 



2t"/xcos ( k 



27T$ 



(5) 



where k e ,k h = 2n£/M are the single particle wave vectors in the absence of a magnetic field 



12j, and K = k e + k h . In the continuous limit, when M — > 00, d — > but the radius of 
the ring is fixed, and considering the energy of the ground state in the absence of the field 
for the electron and hole as zero, Eq. (JHJ) reduces to the standard parabolic expressions for 
charges in a continuous ring threaded by a magnetic flux and no interaction 



B. Exact results 



For 7 7^ 0, the ground state in each momentum sector associated with (J3J) is given by 




Q 

(6) 
(7) 
(8) 
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where 

me 

c n = si ^ [Mu] {sinh[(M - n)v) + e^ Me sinh(^)} , (9) 

9 = arg(g), cosh(^) = —E/2\q\ and cq a normalization factor. Details of the calculation can 
be found in appendix I A II The energy spectrum E v = —2\q\ cosh(z/) is determined by the 
solutions of the equation 



tanh(Mz/) 2\q\ 



sinh(i') 7 



cos(Mfl) 
cosh(Mz/) _ 



(10) 



As we show in appendix IA 21 bound state solutions to (fTU|) exist iff 



7>^(l-cosM#). (11) 

These solutions correspond to the bound electron-hole pair, namely the exciton. In the limit 
M — > oo, d constant, the energy of the exciton is = — a/7 2 + 4|g| 2 , independent of 
the value of the magnetic field. This value corresponds to the exciton binding energy in a 
straight wire ^. Also, the anharmonic/interaction parameter 7 must scale as 1/d in order 
to support a bound state solution; it goes to infinity as expected for a point-like 5-function 
interaction 



C. Approximations for a large ring 

If Md is large, we can derive some explicit approximate expressions (see appendix IA 31 
for details). First, the energy of the bound state is given by 

2 7 2 



E = E n 



— Marcsinhf -3— 



1 • —coMM0)( """" ; ~ 
Et, 



(12) 



This formula is close to the perturbative exciton energy obtained for a large ring radius in 
the continuum model Also, if d is small enough, condition (jllj) can be approximated as 
2%r(3 > 1 — cos(27r$/$o)> where (3 = arcsinh(7/2|g|)/<i is the inverse decay length of energy 
oscillations given by Eq. ()12j) . 

D. Numerical results 



In Fig.Elwe plot the energy of the ground state for several values of the magnetic flux and 
the circumference p = 2nr of the ring. It is easy to see that the magnitude of the AB effect 
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decreases rapidly with increasing ring radius. As expected the maximum amplitude of 
the AB oscillations corresponds to integer or half-integer values of $/$o- If we analyze the 
(unbound) excited states, the contributions of components corresponding to states where 
fermions are localized in different sites of the lattice are significant, and the AB contribution 
no longer decrease exponentially with the circumference of the ring. 




FIG. 3: The Aharonov-Bohm oscillations of exciton energy (in units of £00) corresponding to three 
values of the (single) ring circumference p (in units of d) obtained by solving Eq. IjlOJI numerically 
for 7 = 2, d = 1, e = 1, 11 = 0.2, and k = 0. (o) M = 5 (p = 5.34), (□) M = 7 (p = 7.24) and (*) 
M = 9 (p = 9.19). The lines denote the approximate solutions given by Eq. (|12|) , 



For the ground state, if we plot the amplitude of the AB oscillations AE = E(ip = 
0)—E((p = 1/4) as a function of circumference p as shown in Fig. HJ we observe an exponential 
decay with a decay length d c which can be calculated numerically by solving Eq. (|10j) . The 
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FIG. 4: Energy oscillations amplitude of AB oscillations (in units of and on a semi-logarithm 
scale) as a function of the circumference of the (single) ring p (in units of circumferential lattice 
spacing d), obtained by solving Eq. IjlOJ) numerically for 7 = 1, corresponding to a weakly bounded 
state (o), 7 = 2.5, a typical bound state (*), and 7 = 5, corresponding to a strongly bound state 
(□). e = 1, /x = 0.2, and = 0. The continuous lines are approximate solutions given by Eq. (|12j) . 
All logarithms are natural logarithms. 

characteristic decay length is a function of the anharmonic parameter 7 and, as shown in 
Fig. it decreases with increasing 7. 



10 



1 2 3 4 5 6 7 

y 

FIG. 5: Decay length d c of AB oscillations (in units of the circumferential lattice spacing d) for a 
single ring as a function of the anharmonic parameter 7. e = 1, jj, = 0.2, and k = 0, obtained by 
numerical solution of Eq. H1U|) (*). The continuous line is the approximate solution given by Eq. 

ci- 
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IV. THE CASE OF TWO RINGS 



A. The Hamiltonian matrix 



In the case of two rings (N = 2) with M sites on each, it is again possible to find the 
Hamiltonian matrix H^ K ' in each total momentum K subspace, 



H. 



(K) 





H\2 


H13 


Hu 


TT* 

H 12 


H22 


H23 


H13 


TT* 

H 13 


TT* 

H 23 


H33 


H\2 


TT* 

n u 


TT* 

H 13 


TT* 

H 12 


H^A 



(13) 



where each Hij submatrix is of dimension M x M. The diagonal blocks are given by Hn = 
^(7, qi), H u = ^(7, q 2 ), H 22 = .F(0, q 12 ), H 33 = J"(0, q 2 \) with the M x M matrix F(p, q) 
given by 

p q* . . q 
q q* . 
q q* . . 



?{p, q) 



(14) 



. . . q q* 
_ q* . . . q 

The off-diagonal blocks are if 12 = fit^t, H i3 = t^t, where 1 is a M x M unit matrix, and 
the parameters are 



qi = t\e 2m ^' M ^ + e- lK ), 

q 2 = t l 2 l e 2m ^ /M (n + e- lK ), 

q 12 = t\e 2m ^' M e~ lK + t\e 2 ^ M ^ 

q 21 = t\e 2m ^ /M e- lK + t\e 2m ^' M ^. 



(15) 
(16) 
(17) 
(18) 



Other submatrices are null. Clearly, if t± — and r\ = r 2 , the system reduces to the ID 
case studied before. 

If Ar is small enough, the dominant elements of (jlHjl are H12 and H\ 3 and standard 
perturbation theory can be applied. In first order approximation, the spectrum of H 2 
is given by the superimposition of four sub-spectra. These sub-spectra are given by the 
spectrum of —T{2^, Q)/4 shifted by four different values: (1 -fyu)^, (1 — /J,)t^, (/i — l)t^ and 
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(—1 — , where Q = q% + q 2 + q±2 + 921- According to this result, we expect four bound 
states to exist but some of them can be hidden within the continuum bands, as shown in Fig. 
El where just two of them appear. Numerical tests confirm that this first order perturbation 
result provides the right width for the continuum bands, but fails to precisely locate the 
bound states outside these bands. It also explains the intra-band pattern observable in Fig. 
El in terms of overlapping bands of contiguous sub-spectra. 

In the general case, it has not been possible to find an analytical expression for the 
eigenvectors and eigenvalues of ()13|). but they can be calculated numerically for different 
values of the control parameters and even large values of the number of sites M. In a similar 
way as in the previous section, we define Eoo as the energy of the ground state when M — > oo 
and B = 0. This quantity is a function of the distance between both rings and its value can 
be determined numerically. 

B. Fast and slow AB oscillations 

In general, for strong attraction 7 > 1, as shown in Fig. the ground state is mainly a 
bound state, where fermions are located together in the same location and on the same ring. 
But in contrast to the ID ring, we also find a contribution of components corresponding 
to states where fermions are in contiguous sites and on different rings as shown in Fig. 
El Thus, the variation of the ground state energy with $ / $0 consists of a combination of 
periodic functions with different frequencies. We find that the main contribution to the 
AB oscillations corresponds to the components where the fermions are located around the 
inner ring, with period $i/<lv There exist also a further oscillation corresponding to the 
contribution of the components where fermions are located around the outer ring, with period 
$2/$i- Due to the difference between hop pin g coefficients, a localization effect around the 
inner ring takes place, and, as shown in pl\ . this effect is enhanced by the nonlinearity. 
Thus, only when the distance between rings is very small, the contribution of the periodic 
term with period $2/^1 is (numerically) observable. In all cases, as in a ID ring, the 
contribution of these terms to the energy oscillations with the magnetic flux decreases with 
the circumference of the ring. 

There exists also a significant oscillation term coming from the contribution of non- 
localized states, where the electron and the hole are in contiguous sites but on different rings. 
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FIG. 6: Main components of the wave function corresponding to the ground state as a function of 
the distance between rings Ar (in units of inner ring circumferential lattice spacing d\). 7 = 2.5, 
e = 1, \i = 0.2, M = 15 and N = 2. ( — ) two fermions in a bound state on the inner ring, (- -) two 
fermions in a bound state on the outer ring, (•— ) two fermions in a non-bound state in contiguous 
sites and an different rings, (*) two fermions in a non-bound state in contiguous sites and on the 
inner ring, (• • •) two fermions in a non-bound state in contiguous sites and on the outer ring. 

These oscillations persist when the circumference of the ring increases, although their period 
increases too. We will refer to this contribution as long-period oscillations. Essentially, the 
main period of these oscillations is given by the phase difference between the magnetic field 
term contribution in each ring. Assuming that 7 = and using Eq. (jSJ), this period is 
given by T = 2M/($ 2 /$i - 1) = 2M/(r 2 /n - 1) = 2Mn/Ar. The amplitude of the long- 
period oscillations decreases when the distance between the rings decreases. In cases where 
these oscillations are dominant, corresponding to cases where the distance between rings is 
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small (or equivalently there is a large hopping coefficients t between rings), the oscillations 
coming from components located around the same ring can be seen as a small modulation 
that decreases with the circumference of the ring. We note that a model in which the hole 
is lo c at e d more towards the centre of the r lng has hee n put forward Q irr the context of 

n 

the self-assembled nano-rings [1]. 

In general, if the distance between the rings is small, the ground state gives AB oscillations 
with period $i/<&o. The first excited state, where the main components corresponding to 
the exciton are localized around the outer ring, oscillates with period $2/^0- If the distance 
between the rings is decreased, the hopping probability between different rings increases 
and the bands overlap. Then an additional oscillation of the energy appears due to the 
non-localized components in different rings. 



C. Flux dependence 

In Fig. we show the flux dependence of the ground state energy for three different 
values of distance Ar and two different values of the circumference of the ring. In the top 
figure we can see the small AB oscillations, corresponding to the existence of the exciton 
state localized mainly around the inner ring. If the circumference is small enough, the AB 
oscillations are thus similar to the ID case. If the circumference is larger, the amplitude of the 
AB oscillations decreases and the long-period oscillations become dominant. If the distance 
between the rings is similar to the distance between neighbouring sites and the circumference 
of the rings is still small enough as shown in the central panel of Fig. [?l the excitonic AB 
oscillations are completely mixed-in with the the oscillations of similar frequencies coming 
from non-bound components in different rings. The bottom figure corresponds to the case 
where the two rings are close to each other. If the circumference is small enough, the 
excitonic AB oscillations appear as a slight modulation of the non-excitonic oscillation due 
to the non-bound components in different rings. If the circumference is large enough, only 
the slow oscillations can be seen. 
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FIG. 7: Ground state energy oscillations of the exciton energy (in units of E function 
of ^i/^o f° r several values of the distance between different rings Ar = r<i — r\ (in units of the 
circumferential lattice spacing d\) and different values of the circumference of the inner ring p\. 
7 = 2.5, e = 1, ii = 0.2 and iV = 2. Ar = 1.5 (top), Ar = 0.7 (centre) and Ar = 0.1 (bottom). 
Continuous lines correspond to M = 5 {p\ ~ 5.34) and dashed lines to M = 15 (p\ ~ 15.11). 



D. Fourier analysis of AB oscillations 

If we analyze the Fourier spectrum of the oscillations of ground state energy as function 
of magnetic flux in cases where the AB oscillations are not negligible (e.g., small M), and 
different values of distance between rings Ar, as shown in Fig. |H1 we clearly detect the 
frequencies corresponding to AB oscillations in the inner ring and the oscillations due to 
non-localized components in different rings. A frequency corresponding to the contribution 
of components around the outer ring can only be seen if the distance between rings is very 
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FIG. 8: Fourier spectrum of the ground state energy oscillations of the exciton energy (in units 
of Eqq) as a function of magnetic flux $i/<I>o f° r = 0.7 (in units of the reticular distance 
di) and M = 5 (pi ~ 5.34). 7 = 2.5, e = 1, \x = 0.2 and N = 2. Note that resonances at the 
frequencies uj\ = 2tt, corresponding to AB oscillations in the inner ring, and u>2 = 2.9, corresponding 
to oscillations due to non-bound components in different rings (as well as their harmonics). All 
logarithms are common (base 10) logarithms. 



E. Increasing the ring separation 

If we study the decay of the amplitude AE of the ground state energy oscillations as 
a function of the distance between different rings, as shown in Fig. El an approximately 
exponential decay due to the AB effect can be detected, when the hopping between different 
rings and the diameter of the ring are small. In this case, the typical decay length is similar to 
the ID equivalent case. The most interesting situation occurs when the hopping coefficient 
between different rings is similar to the hopping coefficient along the inner ring. If the 
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circumference of the ring is small enough, AB oscillations are significant and they decay 
with the radius of the ring. Additionally, long-time oscillations appear whose amplitude 
increases slightly with the circumference (and thus the enclosed flux) of the ring. This leads 
to the non-monotonicity observed in Fig. El 
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FIG. 9: Amplitude of the ground state energy oscillations (in units of E^) as a function of the 
circumference of the ring p (in units of circumferential lattice spacing d), corresponding to different 
values of distance between rings Ar = 0.1 (solid line), Ar = 0.7 (dashed line) and Ar = 1.5 (dot- 
dashed line). 7 = 2.5 e = 1, fi = 0.2 and N = 2. The arrow indicates the region which has been 
shown in more detail in the inset. 



V. GENERAL CASE OF N ANNULAR RINGS 

For arbitrary values of M and N, we have not been able to find a general expression for the 
Hamiltonian matrix. Nevertheless, if M and N are small enough, it is possible to calculate 
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the eigenvalues and eigenvectors by using algebraic manipulation methods to construct an 
exact Hamiltonian matrix in algebraic form and then use a numerical eigenvalue solver for 
each parameter value. 

As shown in Fig. we have found essentially the same phenomena as in the previous 
section. If the distance between the rings is large enough, we can find AB oscillations due 
to the existence of components corresponding to non-bound states where the two fermions 
are in the inner ring. These oscillations decrease very rapidly when the circumference of the 
rings increases. When the distance between rings is small enough, there are also significant 
long-period oscillations due to the contribution of states corresponding to electrons and holes 
located in neighbouring sites on different rings. 

If we study the amplitude of the ground state energy oscillations as function of the number 
of rings N as shown in Fig. El we can observe that there exist a weak dependence with the 
number of rings, but this dependence becomes negligible when N is more than a certain N c , 
which in turn is a function of Ar. If Ar is large enough, N c ~ 1. This suggests that, except 
when the distance between rings is very small, the dependence of the energy of the ground 
state with magnetic field is essentially a ID phenomenon. 

VI. A SYSTEM OF STACKED RINGS 

In the case of a stacked array of rings with identical radius located along the magnetic 
field direction — and thus all of them threaded by the same magnetic flux — our results are 
similar to the single ring case. Only the oscillations of period T = $i/$o are observed. If 
we further analyze the decay of the AE oscillations with the circumference of the rings, we 
find that it is nearly exponentially and the decay length can be studied. As expected, if Ar 
is large, the stack behaviour is similar to isolated rings. But if the rings are close, the decay 
length increases as shown in Fig. El This suggests that the bound exciton states use the 
additional rings as further channels along which to tunnel when prompted by the increasing 
magnetic flux. 
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FIG. 10: Aharonov-Bohm oscillations of the exciton energy (in units of E^) as a function of 
$i/$o f° r N = 5 rings and several values of the distance between different rings Ar = r^+i — 
and different values of the circumference of the inner ring p±. 7 = 2. 5, di = l,e = l,// = 0.2,fe = 0. 
Ar = 1.5di (top) and Ar = 0.1<ii (bottom). Continuous lines correspond to M = 5 (p\ = 5.34) 
and dashed lines to M = 13 (pi = 13.12). 



VII. CONCLUSIONS 

In this paper we have studied the existence of AB oscillations for an exciton in a class 
of 2D annular models described by an attractive fermionic Hubbard model. In all models, 
when the ring separation becomes small, the ID continuum results can be recovered in a 
limit, and the structure of exact and approximate equations for the AB oscillations in the 
discrete lattice model is similar to those for the continuum model. 

Our results for the vertically-stacked ring system show that the addition of further trans- 
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0.01 





FIG. 11: Oscillations of the energy of the ground state as function of the number of rings N. 
A = (E mx — E m i n )/(2E), where E is the exciton energy for B = 0. / = 5, corresponding to 
pi = 5.34 (top) and / = 11, corresponding to p\ = 11.15 (bottom). (*) Ar = 0.1, (□) Ar = 0.7, 
(*) Ar = 1.0 and (o) Ar = 1.5. 7 = 2.5, d x = 1, e = 1 and n = 0.2. 



port channels leads to an increase of the decay length for the AB oscillations in the ground 
state. This effect indicates a route to enhancing the possibility of experimentally verifying 
the excitonic AB effect. 

On the other hand, when additional transport channels are added in-plane to give a disc- 
shaped, multi-ring structure as studied in Sections IlII Al W \ the excitonic AB oscillations 
are quickly lost in the oscillations which are due to the various additional flux periodicities in 
each ring. Already for the 2-ring case, our results show quite complicated oscillation struc- 
tures. Nevertheless, even in the general case studied in Section remnants of the original 
AB oscillations survive. These can be seen best when analyzing the Fourier spectrum. 
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FIG. 12: Decay length as function of the distance between rings corresponding to two stacked 
rings (continuous line) for 7 = 2.5 and reticular distance d = 1. The dotted line corresponds to 
the ID ring and has been calculated using Eq. 1121 



Acknowledgments 

The authors are grateful for partial support under the LOCNET EU network HPRN-CT- 
1999-00163. F. Palmero thanks Heriot-Watt University for hospitality, and the Secretarfa 
de Estado de Education y Universidades (Spain) for financial support. We are most grateful 
to our colleague Oliver Penrose for valuable suggestions and discussions during the early 
stages of this research 



22 



APPENDIX A: EIGENVALUES AND EIGENVECTORS IN A ID RING 



We look for eigenvalues and eigenvectors of the matrix given by Eq. (J3J) . The correspond- 
ing eigenvalue equation is 

£|V> = £|V), (Ai) 

where \ip) = (ip , ipi, • • • ,ipM-i) T = J2 m =o ^ml 771 ) w ith suitable (i.e. normalized, complete 
and orthogonal) atomic wave functions at each site m. Due to the special circulant form of 
P| (up to the modified corner element 7) we may consider \ip) as being M-periodic, that is 



1. Eigenvectors and eigenvalues 

The system of equations generated by (|A1|) is 

7^o + q*^i + <pI>m-i = -E4>q, (A2) 

and 

# m _i + q*ip m +l = ~E ip m , (A3) 

with m G {1, • • • , M — 1} with periodic boundary conditions ip m +M = ipm- in addition, we 
choose the phase of ipo to be zero, that is ipo £ ^- This fixes the overall phase factor up 
to which an eigenvector can be determined. We now solve the difference equation (jA3|) . Its 
characteristic equation reads 

q*r 2 + Er + q = , (A4) 

whose solutions are 

r± = (e ± V^l) e ie = e i{e±v) , (A5) 

where we have defined 

E 

q = \q\e l9 ; e = — — — r := cosz/. (A6) 
Note that if |e| > 1, v becomes purely imaginary. We can now write the solution of ()A3j) as 

ip m = [A cos(uM) + B sm(um)] e im6 . (A7) 

It is obvious that A = ipo in the equation above, and using the periodicity condition ipo = ipM 

B = f°. \e~ lM9 - cos(z/M)l . (A8) 
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Thus 

^ = -7\n{e ,mfl sin[(M- m ),] 
sm(z/M) k 

+ e- i{M ~ m)0 sm(um)} . (A9) 

With the convention that ipo G R, we obtain from the latter equation that ipM-m = V'm- 
The normalisation factor f/> can be calculated exactly and reads 

= (A10) 

sin Mv 

r M _ smA^cosM^ + cog m (e^_McosMv)l 1/2 ' 

L tan y V tan v 11 

We now turn to the determination of the spectrum. Eigenvectors of (jlj) have to verify the 
"boundary condition" (jA2|) . Using the general expression derived in (jA9|) . this additional 
constraint results in 

tan(Mz/) _ 2|g 
sin z/ 7 

Given that if is Hermitian, its eigenvalues are real. As we have defined E = — 2|g|cosz/, 
we see that Eq. (|A11|) has to be solved for v G R, v G iR or z/ G 7r + zR (the only distinct 
possibilities for to be real). The solutions on the real axis [y G R) represent extended 
states whereas the unique solution on the imaginary axis [y G iM. or v G 7r + iR ) if it exists, 
represents the bound state. Once Eq. (jAllj) has been solved for z/, the spectrum is obtained 
from E = —2\q\ cosza 



cos(Mfl) 
cos(Mz/) 



(All) 



2. Existence of a bound state 



We restrict hereafter to positive values of 7. It is obvious that the case 7 < can be 
handled from the latter as H(—j,q) = —H( , ~f,e t7T q). For 7 > the eigenvector equation is 
obtained from (jAll|) with v = if], 7] G R. Then, (|A11|) reads 

tanh(Mr/) _ 2\q 
sinh 7] 7 

A necessary and sufficient condition for a bound state to exist is given by Eq. (fTTj) . To prove 
the result we first rewrite ()A12|) as fiji) = 2\q\/j where 

„. . sinhM77 

jvn) = . 

sinh r] (cosh Mr] — cos M9) 
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1 - 



cos(Mfl) 
cosh(M?7) 



(A12) 



Differentiating / with respect to 77 we can show that, for M > 1 and 77 > 0, f'{rj) < 0. 

Now, lim^Q f(j]) = M/(l — cos Mrj) and lirn^oo f(rj) = 0. Therefore, the solution of 
f(rj) = 2\q\/ , ~f exists if 2|g|/7 < M/(l — cos Mr/) which establishes the result. 

Note that when 7 = 2|g|/M(l - cosM6»), the solution of (jMj) is 77 = 0. Then v = 
and ()A9jl shows that the bound state is not exponentially localized because ip m oc e m6 (M — 
m + me~ lMe ) . This limiting case cannot really be considered a bound state, hence the strict 
inequality in (fTTj). 



3. Approximation for large (but finite) number of sites M 

In this case it is possible to derive explicit approximate expressions for the eigenvalues 
both for extended and bound states. 



a. Extended states (v £ M) 

Let us rewrite (jAll|) in the following form 

sin[Mi/ - cp(u)) = - sin <p(v) cos(M0), (A13) 

where tan<£>(V) = 2\q\siai//j. We can restrict our investigation to the interval v £ [0,7r) 
within which we will find M — 1 eigenvalues or M according to whether the condition (jlljl 
for the existence of a bound state is satisfied or not. If M is large, v 9 ( z/ ) varies slowly as 
compared to Mv and the zeroth order approximate solution for (|A13|) is 

^ (0) = ^, /£{0,---,M-l}. (A14) 

Let us seek an approximate solution of the form v\ = + ^ . Then the first correction is 
given by 

= (A15) 
{"P (^ (0) ) + arcsin sin( ^ (^ (0) ) cos(M^) I . 



b. Bound state (v E iR) 

In this case, let v = irj. Then (jAll|) is given by (|1U|). As M — > 00, a first solution is 

7 



r7 (0) = arcsinh ( ^- ) , (A16) 
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and the corresponding energy is 

E {0) = -2\q\ cosz/ = -2\q\ coshr/ 0) = -y/A\q\ 2 + y 2 , 



(A17) 



the solution for an infinite lattice. Notice that this value no longer depends on 9 (i.e. the 
magnetic field). We can now use this last value to solve (|A12|) by a first Newton iteration 
for instance 



V 



where 77W 



fin 



(0)\ 



with 



/(*) 



tanh(Mx) 2\q\ 



sinh x 



7 



1 - 



cos(M#) 
cosh(Ma;) 



As 77 W ~ -2tanhr/°) cos(M6»)e 



we finally obtain 



i£ = — 2|g| cosh 77 
~ -2|g|coshr/ (0) [l + 2tanh2r/ {0) cos(M^) e~ M " (0) 

= - v / 4|g| 2 + 7 2 [i+ 

2 " cos(M^) e - Marcsinh (*i 



(A18) 



(A19) 



(A20) 



4|g| 2 + 7 2 

The factor containing the magnetic field through cos(M8) decreases exponentially with the 
number of sites M, that is with the radius of the ring (the lattice spacing is assumed to be 
constant). 
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